Ground State Energy of Mean-field Model of 
Interacting Bosons in Bernoulli Potential 

Michael Bishop, Jan Wehr 
January 1, 2013 

Abstract 

This paper explores a system of interacting 'soft core' bosons in the 
Gross-Pitaevskii mean-field approximation in a random Bernoulli poten- 
tial. First, a condition for derealization of the ground state wave func- 
tion is proved which depends on the number of particles and interaction 
strength. Using this condition, asymptotics for ground state energy per 
particle are derived in the large system limit for small values of the cou- 
pling constant. Our methods directly describe the shape of the ground 
state in a given realization of the random potential. 

1 Introduction 

In his seminal work [4], Anderson discovered that quantum mechanics behave 
very differently in disordered environments than in periodic environments. Since 
this discovery, there has been a lot of effort both in physics and in mathematics 
to understand this behavior more completely. See [1] which contains a recent 
physical survey of Anderson localization in a context relevant for the present 
paper and [T7] for an excellent introduction to rigorous mathematical results. 

The experimental realization of cold atoms [TTJ[T2] has added to the interest 
in the topic and led to new questions about the role of interactions in quantum 
systems. Physicists have researched the relation between disorder and interac- 
tion [3,6, 13, 15,22]. There is still much debate on the nature of phase transitions 
in both random 24] and quasiperiodic potentials |21j . Systems similar to the one 
discussed in this paper have been experimentally realized |S]. For infinite local 
'hard core' interactions, see [2]. Lieb, Seiringer, and Yngvason have rigorously 
shown convergence of the ground state of a bosonic system with nonrandom 
harmonic potential to the mean- field approximation [10 , 18T-E0] . The dynamics 
of the nonlinear Schrddinger equation with a random potential is studied in [T3] ; 
see references therein. 

Recently, Seiringer, Yngvason, and Zagrebnov considered Bose-Einstein Con- 
densation in systems with randomly placed point scatterers in the continuous 
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setting; they demonstrated the existence of a condensate in such systems (under 
certain conditions on the interaction) as well as a presence of a phase transi- 
tion. [23J . Their analysis is based on a detailed study of the statistics of random 
potential realizations, similar to what we do here. Our model can be viewed 
as a discrete version of the model considered in [53], with the space variable 
rescaled by L and the parameters in the Hamiltonian and the random potential 
scaling as follows: 7 = jqL 2 , a = <jqL, v = vqL. 

For interacting systems in random potentials, Anderson localization com- 
petes with derealization caused by repulsive interaction. It is not obvious which 
of the two mechanisms dominates the system's behavior, even in the ground 
state. One approach to this question is to start from a noninteracting system 
and treat interaction as a perturbation. A system of bosons with no interac- 
tion places each particle in the single particle ground state. When interactions 
are added to a finite- volume system, they can be controlled for small values of g. 

The goal of this paper is to study effects of disorder in a system of bosons 
which interact with a weak repulsive 'soft core' force as the number of parti- 
cles and system size are taken to infinity, proving two results. Theorem 1 is a 
general statement about derealization effects of such interaction, applying to 
both discrete and continuous versions of the model. It is then applied, together 
with detailed analysis of the energy functional, to the one-dimensional system 
with a Bernoulli-distributed potential — Theorems 2 and 3. We derive there an 
asymptotic formula for the ground state energy per particle as the product of 
the interaction strength and the particle density goes to zero. We directly study 
the way the geometry of each realization of the random potential determines 
the ground state wave function. The methods are inspired by an adaptation of 
the technique used in [7J. 

In the discrete setting, the state of N + 1 bosons in the lattice cube A = 
{0, . . . , L+l} d with length L is described by a normalized wave function |$(xi, . . . , xjv+i)), 
a function in the symmetric subspace of <g> N+1 L 2 (A, //) with Dirichlet boundary 
conditions, where the measure fi is the counting measure and X{ are positions 
of the particles in A. The Hamiltonian of the system is given by 



JV+l 




(1) 



where Hi is the single particle Hamiltonian acting on the i-th particle and 
U is the potential of the interaction between particles. The single particle 
Hamiltonian is the (random) Schrodinger operator 



H 



A + V 



(2) 



where A is the discrete Laplacian and V is a (random) multiplication operator, 
which in this paper will always be bounded below (without loss of generality by 
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0). The interaction U is a 'soft core' interaction of the form 



(3) 

where the 5 is the Kronecker delta and the coupling constant g is positive, mak- 
ing the interaction repulsive. This repulsion makes it energetically unfavorable 
for bosons to occupy the same space, but does not exclude the possibility, un- 
like in the the case of 'hard core' interactions where g = oo. Because of the 
difficulty of working in a large tensor space, the multi-particle bosonic states 
are approximated by the Gross-Pitaevskii mean- field wave function [5]. In this 
approximation, each boson is assumed to be in the same state <f> which defines 
the state of the whole system: |$) = \<f>(xi)) . . . \<P(xn+i))- The Hamiltonian 
becomes 

(« + l)H, + 9(A ' + 2 1)(A V . (4) 

a nonlinear random Schrodinger operator (the second term acts as a potential, 
which depends on the state itself — hence the nonlinearity) . The approximation 
exchanges the difficulties arising from Bose statistics for the nonlinearity of the 
new problem. The associated per particle energy functional is 

e w -EE +y)- ^)i 2 + v(*M*)\* + ^l#*)l 4 (5) 

x£A\y\=l 

Remark: In the infinite volume limit, if the particle number increases pro- 
portionally to the system size, the total energy diverges and the natural quantity 
is energy per particle. In this paper, any mention of energy refers to the per 
particle energy. 

In the continuous setting, the state of N + 1 bosons in the box = [0, L] d 
with linear size L is described by the wave function |$(a;i, . . . , xn+i)}, a function 
in the symmetric subspace of (g) Ar+1 L 2 (17, /x) with Dirichlet boundary conditions, 
where the measure /i is the Lebesgue measure, and Xi are the positions of the 
particles in f2. The Hamiltonian of the system is given by 

N+l 

H=Y,H i + y £u(x i ,x J ) (6) 

i=l vjtj 

where Hi is the single particle Hamiltonian acting on the i-th particle and U is 
the interaction potential between particles. The single particle Hamiltonian is 
given by the (random) Schrodinger operator 

H = —A + Vi (7) 

where A is the Laplacian and V is a (random) multiplication operator. The 
potential is bounded below (without loss of generality by 0). The interaction U 
is a 'soft core' interaction of the form 

gS(xi-Xj) (8) 
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where the S is the Dirac delta and g > corresponds to a repulsive interaction. 
In the mean-field aprroximation, each boson is assumed to be in the same state 
<p which defines a state of the whole system: |$) = \<p(xi)) . . . \4>{xn+i))- The 
Hamiltonian becomes 

(N + w + 9(N+ ; m iof, w 

a nonlinear Schrodinger operator. The associated per particle energy functional 
is 



E[$] = ( dx 
Jn 



|^(x)| 2 + ^)l^)| 2 + ^l^)| 4 



(10) 



In section 2, we address the localization and derealization of states in this 
system for both the discrete and continuous settings. For any state cj>, define 
the set 

A >£ (</>) = {x: \<l>{x)\ > J^} (11) 

as the set of points where the absolute value of <f>(x) is greater than its average 
magnitude per site, multiplied by a constant e. This is a natural set in the study 
of localization of low energy states. 

Theorem 1: Assume that V > 0, for both the discrete and continuous 
settings, for a state <j) with energy E' = E[<p], the measure of the set X >e obeys 
the lower bound 

aN(l - e 2 ) 2 

M*>e(0) > 9 \ E , ' (12) 

This theorem implies that Anderson-type localization for low energy states, 
E' w 0, requires strong conditions on gN. For a large number of particles or 
strong interaction, the wave function must fill a significant amount of space, 
meaning the repulsion dominates the Anderson localization effects caused by 
the random potential V. The repulsive interaction, though local, forces overlap 
to become energetically expensive if too many bosons occupy the same place, 
causing the ground state to spread. If a low energy state is localized to some 
length £, gN must be smaller than £ d . To recover an Anderson type localization, 
a necessary condition is gN < 0(E'£ d ). In physical experiments, the interaction 
constant g is a controlled parameter rather than a variable dependent on the 
particle number and the particle density is approximately constant (N w pL d ) 
to ensure the existence of thermodynamic limits. In these systems, a low-energy 
state must occupy a nonzero fraction of the volume rather than being localized 
with some lower order localization length I. 

In section 3, the above theorem is used to describe the ground state of the 
one-dimensional lattice system where the potential V is Bernoulli-distributed, 
the interaction constant g is small, and the particle density is approximately 
constant: N + 1 w pL. The ground state minimizes the energy (per particle) 
given by 

E{4>) = E (V (x)| 2 + V(xMx)\ 2 + ^|#r)| 4 ) (13) 
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where <j> € f 2 {0, . . . , L + 1} with Dirichlet boundary conditions and V is a mul- 
tiplication operator by a function V(x) (of a discrete argument), where V(x) 
are IID Bernoulli random variables: P[V(x) = 0] = p, P[V(x) = b] = 1 — p = q, 
and b > is a constant. 

Definition: The finite volume ground state is defined as the normalized 
state 0q L ' which minimizes the energy functional E[<f>] and the corresponding 
energy, E$ L \ is called the ground state energy. 



Theorem 2: For any g, p, p, and 6, 

lim E ( L) = E (14) 

L— >oo 

where Eq is a nonrandom function of the above parameters. 



After taking the infinite volume limit, we want to understand the behavior 
of the ground state for small gp, when the nonlinear term is taken to zero. 

Theorem 3: For the one-dimensional lattice Gross-Pitaevskii model with 
Bernoulli disorder, the ground state energy E satisfies the following condition 
with probability one. 

< liminf E \ogl(gp) < lim sup E logl(gp) < oo (15) 

Theorem 3 is an illustration of Theorem 1. In Theorem 3, the ground state 
is approximated by sine waves on intervals of zero potential longer than some 
minimum interval length and by zero everywhere else. The ground state energy 
is bounded above by ^ og C ^ gp ^2 ■ Using Theorem 1, p,(X >£ (<j))) is bounded below 

by 

gpL{l-^f{\og p {gp)f 

2C+ 1 j 

which is proportional to the system size. The proof of Theorem 3 shows that 
for the ground state this lower bound is asymptotically accurate. 

Remark: We expect that the upper and lower limit in the above statement 
are equal. One should be able to close the gap between our lower and upper 
bounds using more accurate variational functions — the solution of the discrete 
nonlinear Schrodinger equation on an interval with Dirichlet boundary condi- 
tions, which can be thought of as a discrete version of the Jacobi elliptic sine 
function. 
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2 Proof of Theorem 1 

Theorem 1: If the energy of a state is E' = E[4>] and the potential V is 
nonnegative, then for e € (0, 1) 

,(*>,»> > ^2 (17) 
Proof: The £ 2 -norm of the state <f> restricted to X< e is bounded by 

||^)k<J| 2 <M^<,)^<e 2 

Because the state <f> has L 2 -norm equal to one, 

||^)|x >e ||>l-e 2 

The energy of <f> is bounded below by its interaction energy on X >e . The inter- 
action energy is bounded below using Schwarz's Inequality: 

/ \f\ 2 dS<([ d(,)([ 1/1% 

Jx >e / \Jx >e / \Jx >e 

which bounds the L 4 -norm below by 



L 



lx >e (4>) M^>eW) 

The interaction energy is thus bounded below by 

gN{l-e 2 ) 2 

and bounded above by E' . It follows that 
which gives the desired lower bound: 

^>^l (18) 

and completes the proof. □ 

Remark: By shifting the energy, one can easily generalize the above theorem 
to arbitrary potentials bounded below. If V > V m in , 



2(E> - V min ) 
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3 Ground State Estimates for Weak Interaction 
in Bernoulli Potentials 



Theorems 2 and 3 are specifically for the system on the one-dimensional lattice 
with Bernoulli potential. 

In the proof of Theorem 2, we will apply Kingman's subadditive ergodic the- 
orem, using the version from the standard probability reference [9|; Kingman's 
original formulation |16) would also be sufficient for our purposes. 

Kingman's Subadditive Ergodic Theorem: Suppose X m , n , < m < n 
satisfy: 

(i) Xo, m + X m ^ n > Xo,n- 

(ii) {X nk _( n+1 -) kl n > 1} is a stationary sequence for each k. 

(iii) The distribution of {X m ^ m+ k, k > 1} does not depend on m. 

(iv) E[Xg~-|J < co and for each n, E[Xo, n ] > 7o n i where 70 > —00. 

Then 

(a) limn^oo E[X 0>n ]/n = inf m E[X 0>m ]/m = 7. 

(b) X = lim„_ i . 00 Xo/n exists almost surely and in L 1 , so E[X] = 7. 

(c) If all the stationary sequences in (ii) are ergodic then X — 7 almost surely. 

Proof of Theorem 2: We have 



■0 =1 r— ; ^ 



v =- 

i=0 



where ip € £ 2 {0, ■ ■ ■ L} with Dirichlet boundary conditions. With <j) = yLip, we 



can rewrite E Q as 



L-1 



^ = 7 min V 100" + 1) - 0(,)| 2 + y(j)|0(j)| 2 + ( 2 °) 
which is equal to j^Xq^l, where 

Xo.L= min J2\Hj + ^-Hj)\ 2 + V(j)\m 2 + i ^^ (21) 

11011=^^ 2 

The process Xo,£ satisfies the assumptions of Kingman's theorem. For (i), Xq.l 
is subadditive since the only difference between Xq,m + Xm.l and Xq^l is the 
restriction <f)(M) = in the definition of the former. Properties (ii) and (iii) 
hold because the V(j) are independent. For (iv), X^ s < 2 + 6+^ and X . n > 0. 
Therefore, 

lim E ( Q L) = lim ^ = E Q (22) 

L— foo L— ^-oo Xj 
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almost surely. □ 
Theorem 3 says, in essence, that the limit of as L —> oo to leading order 



C 



ia (i g ( gp ))2 f° r some constant C when gp small. The proof of Theorem 3 does 
give the relationship for small gp, however it does not provide the precise value 
of C . However, the proof of Theorem 3 does provide a good approximation of 
the true ground state. 

Theorem 3: For the one-dimensional lattice Gross-Pitaevskii model with 
Bernoulli disorder, the ground state energy E satisfies the following condition 
with probability one. 

< lim inf E \og 2 Jgp) < lim sup E logl(gp) < oo (23) 

The structure of the ground state depends on the distribution of the inter- 
vals of zero potential. These intervals are colloquially referred to as 'lakes' or 
'islands.' A realization of the potential V is determined by alternating intervals 
of zero potential and positive potential. In a sequence of Bernoulli random vari- 
ables, the lengths of intervals of zero potential are independent and distributed 
geometrically: P[Li = x] = qp x . If the system size L is fixed, the intervals of 
zero and positive potential are not independent; they are subject to the condi- 
tion that the sum of their lengths is exactly L. Their number is not constant, 
but it easy to show that it satisfies a law of large numbers. The considered sys- 
tem can by approximated by a system of variable length, in which the number 
of intervals is fixed at the value dictated by the law of large numbers. Such an 
approximation was carried out in detail in [7] by standard probabilistic methods 
and will not be repeated here. 

Let us thus fix the number of independent intervals of zero and positive 
potential to be 2n. This makes L = J27=i(Li + Li) a random variable, the sum 
of n geometrically distributed intervals of zero potential and n geometrically 
distributed intervals of positive potential . There intervals of zero potential 
have lengths Lj with the distribution 

P[U = x] = qp*- 1 (24) 

for integer values of x greater than or equal to one (an interval must have at 
least one site of zero potential). Likewise, lengths of there intervals of positive 
potential are distributed according to 

P[Li = x] = pq x ~ l (25) 

The variable i indexes the intervals and takes values 1, . . . , n. The total system 
size L is the sum of the random interval lengths. The system size L has expected 
value 



- (26) 
pq 
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Remark: By the Law of Large Numbers [S], with probability one in the 
limit n — > oo, the difference of both L and ^2 L >x Li and their expectations 
has order less than n. These controls occur with probability one in the limit 
n — > oo, further referred to as: 



where 



E 



\L--\=o{n) 
pq 

| £ L«-E[£ Li]\ = o(n) 

Li>x Li>x 

J2 u = n J2 y p ^ = y\ 

,Li>x J y> La:J +1 

pq\ J 



(27) 
(28) 



(29) 



The result of Theorem 3 is obtained using the following strategy. An upper 
bound on the ground state energy can be generated by evaluating the energy 
functional on any test function. The method is to find a test function with 
asymptotics similar to a demonstrable lower bound. The process is iterative; 
first, a test function is evaluated to find an upper bound on the ground state 
energy. Second, this upper bound is used on the norm of the true ground state 
restricted to various sets, such as the set of sites of positive potential, in or- 
der to isolate the sites where the £ 2 -norm of the ground state is concentrated. 
Third, the energy of the ground state on these sets is minimized to find a lower 
bound. This lower bound may give intuition for a better choice of test function 
for the upper bound. The process repeats until the upper and lower bounds on 
the ground state energy are asymptotically similar. The proof of Theorem 3 is 
simply the final iteration of this process. 

Proof of Theorem 3, Upper Bound: The energy of any test function 
bounds the ground state energy from above. Consider the test function ip defined 
as follows: for an interval of zero potential with length Li, the function is a sine 
wave rrii^J -^pj sin( ), where nii is the L 2 -norm of the function restricted 

to the interval. On intervals of high potential, ip is zero. For intervals of zero 
potential with length L; > \og p (gp) + log p (log p (gp)), we let 

m? - ^ T (30) 

2^Li>log p (gp)+log p (log p (gp)) 1 

and for shorter intervals, = 0. This makes the kinetic energy term and 
interaction energy term have the same asymptotic order as gp — > 0. This test 
function also satisfies the normalization criterion J^. m| = 1. The kinetic energy 
of the discrete sine wave on a specific interval of zero potential is bounded above 
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by (i^+iy} 0- The interaction energy, J2 X \\(f>(x)\\ 4 is equal to 39 ^ m ' . 
Summing the upper bound on kinetic energy and the interaction energy of test 
function ip over the space, the total energy of the test function ip is bounded 
above by 

SgpL , ^ 



4 EL 1 >iog p ( S p)+iog p (iog p ( SP )) L * (logp(gp) + log p (log p (gp)) + l) 2 

where the second term is an overestimate of the kinetic energy, treating each in- 
terval as the shortest interval admitted. Both L and V r ^, , /, , ,\L; 

depend on the realization of the potential, but by equations (27-29), with prob- 
ability one in the limit n — > oo, 



Ti 

L = — + o(n) (32) 
pq 



and using [x\ > x — 1, 



£;>l°gp(3P)+l°gp(l°gp(9P)) 

> — 5plog p (sp) (p + qlog p (gp) + qlog p (log p (gp))) + o{n) (33) 

by equation (31). The interaction energy is bounded in the limit as n — > oo 
with probability one by 



4glog p (5P) [logp(sp) + log p (logp(gp)) + | 

In the limit as <?p — > 0, the leading order term of the upper bound on the ground 
state energy is 

a 



logp(ffp) 

with C = ^ + 7T 2 . Thus, limsup^^ E logp(gp) < oo. □ 
Proof of Lower Bound: 

For each n and gp, the ground state </>o is well-defined but not explicitly 
known. The lattice will be partitioned into four sets: the set of sites of high 
potential, denoted the set of sites on intervals of zero potential longer than 
log p (gp), denoted Mi ong ; the set of sites on intervals of zero potential shorter 
than log p (gp) where the kinetic energy cannot be easily bounded below, denoted 
M, 

ightl and the set of sites on intervals of zero potential shorter than \og p (gp) 
where the kinetic energy can be easily bounded below, denoted Mheavy 

The lower bound is shown as follows. The ground state energy is bounded 
below by a lower bound of the ground state energy restricted to Mheavy The 
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kinetic energy for a given interval in Mheavy is bounded below by Lemma 1. 
The interaction energy is bounded below by the bound of Theorem 1. Lemma 2 
provides a lower bound on the norm of <po restricted to Mheavy, which converges 
to one in the limit gp — > 0. Using Lagrange multipliers, the lower bound on 
kinetic and interaction energy on each interval is minimized over the m,-, the 
norms of restrictions of the state to each interval. This minimization determines 
a minimal interval length for any interval in Mheavy The number of sites in 
Mheavy is estimated above by the number of sites on intervals longer that this 
minimal interval length. Using this and the lower bound of the norm of 0o re- 
stricted to Mheavy , we obtain the desired lower bound of the ground state energy. 

Without loss of generality, the ground state wave function can be assumed to 
be non-negative. By a standard argument: if the ground state does not have the 
same complex phase at each site, the kinetic energy can by reduced by making 
the phases equal. Since the potential and interaction energy only depend on the 
magnitude of the state at each site and not on the complex phase, the energy 
of the resulting state is strictly smaller. Therefore the ground state must have 
the same complex phase and can be assumed to be positive. 

To define Mu g ht and Mheavy, the ground state energy on an interval of zero 
potential is approximated by that of the the minimizer of the kinetic energy on 
the interval — the discrete sine wave. For a given interval with length Li, the 
ground state determines boundary values rriiSf and wiiS* 1 on the sites of high 
potential to the left and right of the interval, respectively, where rm is the norm 
of the ground state on the interval. The boundary values are positive by the 
above argument. The kinetic-energy-minimizing state is of the form 

^=FT s '° ( I-TT + " ) (34) 

where the sine wave is normalized to m, by Ci € [1, >/2], stretched by Si £ (0, 1), 
and shifted by U; all three are determined by the S^s and rrii [7]. Heavy inter- 
vals have relatively small S^s which determine a lower bound on kinetic energy. 
Light intervals have large Si 's which do not admit a good lower bound on kinetic 
energy, but do allow an upper bound on the norm of the ground state on this 
intervals. 

Definition: For intervals of zero potential with length less than \og p (gp), 
an interval is in Mheavy if for the ground state cf>o, 

max(^,<5f)<-i= (35) 
zy Li 

An interval is in Mu g ht, 

max(^,<5f)>-l= (36) 

l\l Li 

These definitions can be restated using m, and are thus directly determined 
by the values of the ground state wave function on the sites adjacent to the 
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zero potential intervals. The definition is stated above without m, in order to 
separate the shape and curvature of the sine wave from its norm and to facilitate 
estimates. The norm of 4>o restricted to a heavy interval is bounded below by 



m l > m^^/Li max(<5f , 5f ) > m l vZ" max(5f , 5f ) (37) 

where the last term is the norm of the constant function rrii max(5^, 5f~) on the 
interval. This means that a sine wave with this norm achieves its maximum 
rather than being nearly flat. 

Lemma 1: The kinetic energy of the ground state restricted to an interval 
in Mheavy is bounded below by 

'"■('-^(I^TF (38) 

Proof of Lemma 1: The kinetic energy of a heavy interval is bounded 
below by the kinetic energy of the sine wave with norm vtii and boundary con- 
ditions rriiSf and rmSf . The energy of the minimizer J^"^ sm ( l^+i + *») * s 
mf sin 2 ( L Si ^ r 1 ). To solve for Si, note that the function must satisfy the boundary 
conditions 

: sin(ti) = rmSf (39) 



sin(si7r + ti) — rmSf- (40) 



The left boundary condition is solved using the inverse of sine on [0, arcsin(ir) 

I 2 ' 2 



The right boundary condition is solved using the inverse of sine on [~, 4?], 



arcsin(x) + ir. We obtain as in [7]: 



u = arcsin i h ^ u ± 1 1 (4i) 



c, 



Sin + ti = 7T — arcsm I (42) 



Solving for Sj 



1 - I (arcsin (^H) + arcsin ^NEH) ) , 4 , , 



Since arcsin(6») < ^ for 9 < 1, 



Si > 1 - 



>!-J= (44, 
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where c» > 1, Lj > 1, and the definition of a heavy interval is used. □ 
Lemma 2: In the limit as n — > oo, with probability one, 



46 log p (fifp) 



where O(-) is taken with respect to the limit — > 0. 

Proof of Lemma 2: Because the potential energy of the ground state must 
be less than the upper bound on the ground state energy, the norm of the ground 
state restricted to high potential is bounded as follows: 

II4>o1mJ 2 < - ° + f 2 (46) 
H lo g P (.9P)) 

An upper bound on thchc norm of the ground state on intervals longer than 
log p (gp) follows from the fact that the interaction energy must be bounded above 
by the upper bound on the ground state energy. The minimum of the interaction 
energy depends on the number of sites occupied and the norm restricted to the 
set of these sites. It was shown in Theorem 1 that for \\<f>\\ = to', the minimum of 
X)f=i I^WI 4 i s ^TT'i this result is also due to the inequality between arithmetic 
and quadratic means. If to is the norm of 4>o\M long , then the interaction energy 
for the sites in Mi ong is bounded below, using Theorem 1. 

9pLm 4 9pm\f q - o(n)) 

E Li> ^ p{gp) Li- ^( P + pqlo gp (gp))+o(n) 1 j 

Using the upper bound on the ground state energy, the left hand side of the 
above inequality is bounded above by -, — — - — r^. In the limit as n — > oo, 

(log p (sp)j 

/ pqC+ pC+ . . 

" V bg^) + {\og p {gp)Y ( } 



\\<P0\M lo 

For light intervals, 



TO? 



4L _ <m?max(5f,jf) 2 (49) 



The upper bound on the norm of the ground state restricted to sites of high 
potential bounds the norm on the boundary points 

V mj max(#, 5? f < 2 ° + 2 (50) 

b(log p (gp)) 

where the extra factor of 2 is included to cover the cases where two intervals of 
zero potential are separated by a single site of positive potential. The bound 
on the norm of the ground state restricted to light intervals, which by definition 
are shorter than \og p (gp), is 
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E 

light intervals 



light intervals 



!og p (5P) 



hp/ii intervals 



< 



4b (log p (gp)) 

This means that the norm on light intervals is bounded above by 

Sp ( 



\\<h\ 



\M li9ht \\ < mogp{gp) 

The normalization condition requires 



ll^olMjl' + ll^olAf^Jl' + ll^olAf^JI 2 



y 0|M h 



1 



(51) 



(52) 



(53) 



Using the upper bounds on the three other terms, the norm of the ground state 
restricted to Mh eavy gives the desired lower bound 



ll0o|M heovi) 



> 1-0 



y / !°gp(5p) y 



o 



1 



b\og p (gp) 



The energy for a heavy interval is bounded below by 



gpLmi , 2,, 
• mj 1 



2Li 



V27 (£< + i) 2 



□ 



(54) 



where the first term is the minimum of the interaction energy and the second 
term is the minimum of the kinetic energy. Using Lagrange multiplier method, 
the choice of rrii, where i labels the heavy intervals, which minimize this lower 
bound 



E 



gpLmj n ( , 1 \ 2 7r 2 



+ mf 1 _ , , . 



(55) 



under normalization constraint 



E m ? = 1 -° 



yiogp(sp), 



must satisfy 
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dm; 



E 



2L,. 



O 



+ mf 1 - -= -o 



1 

3 



2 2 
7T 



(56) 
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This equation implies that 

nii = (58) 

or 

U ( ( 1 \ 2 7T 2 ^ 



^"SE^-l 1 -^ IfJ (59) 

For an interval to contribute to the minimization of the lower bound, the kinetic 
energy of a heavy interval must be less than A. For the kinetic energy to meet 
this bound, the heavy interval must have length 

It follows from Lemma 2 that the normalization condition requires A to satisfy 
£ i* =1 _ ( ij.of *A ,01) 

Using equation (29) and approximating [a;J , 

V ^=(-+0^)) (WqpW+i+pW+A (62) 

Li>x \P9 1 

which requires 

( ^+o(n)) * apt 1 "*)* +P ( 1 -^)^) « 1 

9p(f q +o{n))\P<l V VV J 

(63) 

The asymptotic behavior of the parameter A is determined by gp. If A is constant 
or taken to infinity, the normalization will not hold because the left side of (63) 
goes to infinity. If A converges to zero, the dominant term is the exponential. 
The correct asymptotic solution for A as gp — > is 



(64) 



This substitution will make left hand side (63) converge to a constant. Therefore, 
the heavy intervals must be longer than log p (<?p) + log p (log p (gp)) + 0(1). The 
size of Mneavy is bounded above by the number of sites on intervals longer than 
this lower bound. This number is bounded above by 

ngp\og p {gp) (log p (0p) + log p (\og p {gp))) + o(n) (65) 
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The interaction energy of any state supported on this number of sites is bounded 
below by 

9P(%+0W) 

^ngp\og p {gp)[\og p (gp) + log p (\og p (gp))} + o(n) 

The lower bound follows. □ 
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